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MAGNETIC  INSULATION  IN  COAXIAL  TRANSMISSION  LINES 


WITH  AN  EXTERNAL  MAGNETIC  FIELD 

Michael  A.  Mostrom,  Michael  E.  Jones,  and  Lester  E.  Thode 

Intense  Particle  Beam  Theory  Group 
Los  Alamos  Scientific  Laboratory 
Los  Alamos,  New  Mexico  87545 

Abstract 

A  relativistic  Brillouin  flow  equilibrium  model  is  assumed  for  the 
plasma  sheath  and  a  cylindrical  analytic  solution  including  a  uniform 
external  magnetic  field  is  obtained,  A  free  parameter  appearing  in  all 
previous  magnetic  insulation  theories  is  shown  to  be  accurately  deter¬ 
mined  by  maximizing  the  transmitted  power  flow.  The  theory  is  supported 
by  particle-in-cell  simulations. 
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The  demands  placed  on  high-voltage  (multi-megavolt)  technology 
have  led  to  a  great  deal  of  interest  in  the  phenomenon  of  magnetic 
insulation;  i.e.,  the  ability  of  magnetic  fields  to  partially  or  com¬ 
pletely  insulate  the  space  between  two  conductors  held  at  a  potential 
difference  far  exceeding  the  threshold  for  field  emission.  Magnetic 

insulation,  as  applied  to  the  case  of  transmission  lines,  has  been 

1-2 

studied  theoretically  by  several  authors  ,  but  in  all  cases  a  param¬ 
eter  (e.g.,  the  radial  extent  of  the  plasma  sheath)  has  been  left 
unspecified.  The  only  exception  has  been  early  particle-in-cell  simu¬ 
lations,  with  no  external  axial  magnetic  field,  which  indicate  that 
magnetic  insulation  operates  at  the  minimum  possible  current  (in  only 
partial  agreement  with  the  results  here).^  Thus,  these  theories  offer 
no  prediction  for  the  actual  operating  current  of  a  magnetically 
insulated  transmission  line.  In  addition,  only  in  the  planar  approx¬ 
imation  has  the  effect  of  an  external  magnetic  field  been  included.^ 
A  complete  theory  of  magnetic  insulation  in  a  transmission  line  is 

necessary  for  the  proper  design  of  experiments  requiring  the  generation 

4  5 

of  intense  relativistic  electron  beams.  ’  It  is  particularly  impor¬ 
tant  that  the  theory  include  the  effects  of  an  external  magnetic  field 
which  is  crucial  to  the  operation  of  foilless  diodes.^  A  foilless 
diode  now  appears  to  be  a  promising  device  for  producing  very  intense, 

highly-collimated,  annular  relativistic  electron  beams. ^  Such  a  beam 

17  20-3 

may  be  capable  of  producing  a  10  -10  cm  ,  multi-kilovolt  plasma 

8 

with  sufficient  power  density  to  implode  a  small  cylindrical  liner. 

In  the  theoretical  analysis  here,  a  complete  theory  of  magnetic 
insulation  in  a  cylindrical  coaxial  transmission  line  is  presented  by 
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using  a  relativistic  Brillouin  flow  model  (zero  radial  velocity)  for 

the  electron  sheath.  This  is  valid  for  the  usual  case  where  the 

voltage  rise  time  is  long  compared  with  the  cyclotron  period  of  an 

9 

electron  in  the  combined  external  and  self -magnetic  fields.  In  any 

event,  the  magnetic  insulation  solution  appears  to  be  insensitive  to 

2 

the  particular  flow  model.  The  theory  presented  here  differs  from 
previous  analysis  in  two  respects:  (1)  a  uniform  external  magnetic 
field  directed  along  the  transmission  line  is  included;  (2)  the  free 
parameter  in  the  previous  theories  is  set  here  by  maximizing  the  power 
flow  while  holding  constant  the  voltage  of  the  wave  propagating  down 
the  transmission  line.  A  very  accurate  cylindrical  analytic  solution 
is  found. 

This  problem  was  also  studied  using  the  two-dimensional,  time- 
dependent,  fully  relativistic  and  electromagnetic,  particle-in-cell 
simulation  code  CCUBE.^^  Simulations  were  carried  out  at  a  variety  of 
voltages,  ratios  of  anode  to  cathode  radii,  and  uniform  external 
magnetic  field  strengths. 

In  the  theoretical  analysis,  an  equilibrium  state  with  no  azi¬ 
muthal  (9)  or  axial  (z)  variation  and  no  radial  velocity  (v^)  is 

assumed.  Then  the  conservation  of  energy  and  canonical  momenta 

2 

immediately  gives  three  integrals  of  the  motion:  e(t)/mc  =  Y"!  >  eA^/mc 
=  YP^,  and  eA^/mc  =  yPqj  where  the  electrostatic  ((f))  and  vector  (A^,Aq) 
potentials  and  electron  velocities  p^  ~  ^  assumed 

to  vanish  at  the  cathode  radius  r  =  a,  and  Y  =  ^ 

and  m  are  the  electron  charge  and  mass,  and  c  is  the  speed  of  light. 
The  only  remaining  quantities  to  determine  are  p^  and  p^.  Amperes  law, 
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Poisson’s  equation,  and  the  above  constants  of  motion  can  be  combined 
to  yield^^ 

=  p2^1)/PY^  (la) 

[InCpy^pg)]  =  Pq/P^Pq 

where  the  prime  denotes  d/dp,  p  =  r/a,  and  p"(l)  =  P"(p  =  1)  is  a 

'Zi  Zt 

constant. 

At  present,  we  know  of  no  exact  analytic  solution  to  Eqs.  (la)  and 
(lb).  The  simulations,  however,  demonstrate  that  p  “r  very 

nearly  independent  of  p.  Taking  r|(P)  “  n(0,  ^  constant,  then  Eq.  (la) 
gives 

y  =  cosh  [ofp^(l)lnp]  (2a) 

=  a  ^  tanh  [ap^(l)lnp]  =  Pg/r)  (2b) 

2  _  2 

where  a  =  1  +  q  .  Equation  (lb)  in  turn  allows  n(P)  to  be  estimated: 

n(p)  =  n(l)exp{0.5[0fp^(l)l  ^  sinh^  [ap^(l)lnp]} 

2 

or,  for  small  p,  n(p)  -  n(l)^^p[  (l^^P)  /2].  Thus,  the  approximation  q  = 
q(l)  is  good  if  the  plasma  sheath  remains  relatively  close  to  the 
cathode  such  as  will  occur  for  large  external  magnetic  fields,  large 
applied  voltage  (as  discussed  later),  or  sufficiently  small  ratio 
(b-a)/a  where  b  is  the  anode  radius.  The  approximation  q  =  n(l)  is 
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also  applicable  for  sufficiently  small  external  magnetic  field,  because 
then  Pq  is  negligibly  small  and  Eq.  (lb)  becomes  irrelevant. 

Using  the  approximation  r)  =  n(l)»  Eq.  (2)  applies  in  the  region 
1  S  p  ^  pp  where  p^  designates  the  outer  radius  of  the  plasma  sheath. 
Outside  the  plasma  sheath  (p^  =  p  =  b/a),  the  electrostatic  potential 
is  given  by  y  =  [yQln(p/Pp)-yp  ln(pa/b) ]/ln(b/aPp) ,  where  =  yCpp) 
and  Yq  ”  Y(P  ”  b/a)  is  the  anode  potential.  Continuity  of  the  radial 
electric  field  at  using  this  expression  for  y  and  Eq.  (2a),  results 
in 

Inpp  =  ln(b/a)[l-yp(yQ-yp)/(yp-l)G] 

G  H  [cosh“^(yp)  +  (yQ-yp)/(yp-i)‘^]  Yp/Cyp-i)'^ 

Also,  Eq.  (2a)  implies  P"(l)  =  cosh  ^(y  )/o(lnp  .  Using  =  -dA  /dr, 

2  P  P  D  Z 

the  total  axial  current  (cathode  boundary  current  plus  electron  sheath 
current)  is  ~  ^TTmc/ep^. 

2 

So  far  the  total  anode  voltage  =  (Yq  “  'l)mc  /e  has  been  treated 
as  a  parameter.  Actually,  in  any  experiment  or  time-dependent  simula¬ 
tion  is  not  a  constant  but  is  a  function  of  the  plasma  sheath  prop¬ 
erties  which  vary  with  time.  This  is  because  the  plasma  alters  the 

transmission  line  impedance  and  induces  reflections.  In  practice  a 

2 

wave  (e.g.,  right-going)  of  known  fixed  voltage  »  (Yj^  *  l)mc  /e  is 

launched  from  a  driver  of  impedance  into  the  transmission  line.  The 

total  steady-state  voltage,  including  reflections,  is  then  related  to 

y  and  the  total  current  in  the  transmission  line  by 
K  U 


(3a) 

(3b) 
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Yq  =  1  +  (Yjj  -  1)(2  -  I)  ,  (A) 

where  I  =  In/I„  and  =  V_/Z_,  so  that  determining  the  steady-state 
UK  K  K  U 

current  I  gives  the  appropriate  In  the  following  discussion,  is 

chosen  to  be  the  vacuum  impedance  of  the  transmission  line, 
(ljQc/27t)£n(b/a) . 

Using  the  above  expression  for  along  with  Eqs.  (3a)  and  (4), 

one  finds 

I  =  (yJ-1)^  G/a(Yp-l)  •  (5) 

p  K 


The  next  step  is  to  determine  a.  For  1  S  p  ^  the  azimuthal 
vector  potential  is  Ag  =  yP^mc/e,  while  for  p^  S  p  ^  b/a,  it  is  Ag  = 


2  2 


(p  -p^)k/p.  The  constants  p^  and  k  are  determined  by  matching  A^^  at  p 


1 


0 


2  2. 


and  by  setting  Ag(p  =  b/a)  =  (b  -a  )Bg/2b  from  conservation  of  axial 
magnetic  flux.  Note  that  the  uniform  external  axial  magnetic  field  Bg 
has  made  its  first  appearance  in  the  theory.  Finally,  matching  the 


axial  magnetic  field  at  p^  from  the  above  two  expressions  for  Ag 
requires  that 


a  =  [l-4<QU^(b/a-a/b)^/(yp-l)]“’^  (6a) 

U  H  X„[l  +  X^+(l-X^)G/ln(b/a)]’^  (6b) 

K  K  K 


where  v  =  P  a/b  and  -  eB^a/mc. 

^p  ^  p  ^0  0 

The  coupled  transcendental  equations 


solution,  with  y^  still  a  free  parameter, 
numerically.  For  Bg  =  0  this  solution 


(3-6)  form  the  basic  analytic 
and  must  in  general  be  solved 
sinijilifies  to  that  of  Ref.  2. 
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The  reason  is  as  yet  unspecified  is  that  only  the  equilibrium 
equations  have  been  solved.  No  indication  of  which  equilibrium  (out  of 
a  possible  continuum)  the  full  time-dependent  equations  would  even¬ 

tually  lead  to  has  yet  been  given.  Whether  this  evolution  through 

various  other  possible  equilibrium  states  results  from  instabilities  or 

12 

statistical  probabilities  is  not  discussed  here,  but  the  time- 

dependent  simulations  performed  here  demonstrate  that  the  system  does 

evolve  to  and  oscillate  about  a  unique  steady-state.  The  oscillations 

could  indicate  that  the  system  is  trying  to  maximize  or  minimize  some 

quantity.  Furthermore,  as  increases  1  decreases  toward  unity  until, 

for  larger  than  some  critical  value,  1  remains  clamped  at  unity. 

This  is  shown  in  Fig.  la,  which  plots  the  operating  impedance  Z/Z^  ^ 

V^/I^Z^  as  a  function  of  for  b/a  =  1.853  and  two  values  of  (1.5 

MeV  and  6.4  MeV).  Thus,  the  system  appears  to  be  maximizing  some 

function  of  I,  P(I),  subject  to  some  constraint  on  the  range  of  1.  As 

Bq  increases,  the  constraint  is  relaxed  until  I  can  reach  closer  and 

closer  to  unity  where  P  attains  its  maximum  value.  The  theory  indeed 

shows  that  I(Yp)  has  a  range  of  values  that  moves  toward  and  eventually 

encompasses  unity  as  B^  increases  from  zero.  A  physical  quantity 

having  the  desired  form  is  P(I)  =  1(2-1),  the  transmitted  power  flow 

2 

down  the  transmission  line  in  units  of  the  launched  power  P^  =  Vj^/Z^. 
It  is,  therefore,  proposed  that  magnetic  insulation  operates  in  the 
maximum  power  flow  mode.^ 

By  this  it  is  not  meant  that  1  =  1,  but  rather  that  is  chosen 
so  that  iC'/p)  is  as  close  to  unity  as  the  constraints  allow.  For 
sufficiently  small  B^,  i(Yp)  ^  ^  ic>r  all  y^  so  that  in  this  regime 
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power  flow  maximization  is  equivalent  to  current  minimization  (but 
remember  is  held  fixed,  not  V^). 

Using  this  argument  to  set  in  Eqs.  (3-6),  the  theoretical  Z/Z^ 
is  plotted  as  a  function  of  in  Fig.  la.  Also,  plotted  in  Figs.  lb,c 
are  the  theoretical  and  simulation  results  for  p  and  the  ratio  I^/I^ 
of  plasma  to  total  current  as  functions  of  Eq.  (2)  predicts  Ip/^Q  “ 
(Vp'"l)/Vp>  which  extends  the  result  of  Reference  (2)  to  include  B^^O. 
Due  to  temporal  and  spatial  variations  and  limits  on  code  resolution, 
the  uncertainties  of  the  simulation  measurements  are  on  the  order  of 
±  5%  for  Z/Z^  and  and  ±  10%  for  Ip/lQ-  In  all  cases,  the  theory  and 
simulation  are  in  close  agreement. 

Several  interesting  observations  can  be  made  from  the  theory: 
(1)  Z/Zq  is  independent  of  Z^  only  for  =  0  in  which  case  Z/Z^  = 
0.82[(y  -  1)/(V  +  1)]^;  (2)  pp  and  Ip/lQ  first  increase  as  B^ 

increases  from  zero  but  eventually  fall  abruptly  as  B^  exceeds  the 
critical  value  where  I  =  1;  (3)  the  plasma  sheath  radius  (maximized 

over  all  B^)  decreases  (i.e.,  the  effective  insulation  increases)  as 
increases;  and  (4)  a  plot  of  the  total  (radially  integrated)  energy  per 
unit  length  along  the  transmission  line  as  a  function  of  does  not 
reveal  any  consistent  pattern  (e.g.,  an  extremum),  valid  for  all  and 

that  could  be  tested  as  an  alternative  principle  for  determining 
the  operating  Observation  (3)  can  be  interpreted  as  due  to  the 

electrons  needing  less  radial  distance  (for  higher  V^^)  to  acquire  the 
necessary  axial  kinetic  energy  and  hence  axial  current  and  required 
to  insulate  themselves.  Observation  (2)  can  he  interpreted  as  due  to 
the  diversion  of  axial  energy  into  azimuthal  energy  making  it  more 
difficult  to  geiH'rate  an  axial  current;  for  large  B^^,  however,  the 
axial  magnetic  field  excc'cds  and  dominatt's  tiu'  insulation. 
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In  conclusion,  a  complete  theory  (i.e.,  no  free  parameters)  of 
magnetic  insulation  in  a  coaxial  transmission  line  with  an  arbitrary 
axial  external  magnetic  field  has  been  presented.  The  analytic  theory 
consists  of  the  coupled  transcendental  equations  (3)’(6)  combined  with 
the  maximum  power  flow  condition  that  |I(Yp)  ’  1 1  be  minimized  over  all 
to  determine  the  actual  operating  value  of  Yp*  Two  crucial  points 
were  necessary  in  deriving  this  result:  (1)  the  launched  voltage  is 
the  proper  quantity  to  hold  constant  (not  the  total  voltage  V^)  when 
one  calculates  for  different  Yp  above;  (2)  including  finite 

external  magnetic  field  effects  allows  the  solution  to  be  studied  as  a 
function  of  an  additional  variable  (^q)  »  comparing  this  variation 

with  simulation  results  leads  to  the  maximum  power  flow  condition. 

The  authors  wish  to  thank  B.  B.  Godfrey  and  M.  C.  Clark  for  their 
helpful  comments. 
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Figure  1 


FIGURE  CAPTION 


Magnetic  field  variation  of  magnetic  insulation  impedance 
Z/Zq,  plasma  sheath  radius  p^,  and  ratio  of  plasma  to  total 
current  I  /!„.  Solid  lines  are  laminar  flow  theory  with 

p  0 

power  maximization.  =  eB^a/mc  =  0.59  (kG)  a  (cm), 


b/a  =  1.853. 
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